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APPLICATION  OF  ADVANCED  FRACTURE  MECHANICS  TECHNOLOGY  TO  ENSURE 
STRUCTURAL  RELIABILITY  IN  CRITICAL  TITANIUM  STRUCTURES 


H.  A.  Ernst  and  J.  D.  Landes 
Materials  Engineering  Department 


ABSTRACT 

This  final  report  describes  the  work  done  in  a  program  designed 
to  assist  the  Navy  in  developing  and  applying  fracture  mechanics 
technology  to  assess  structural  reliability  in  critical  applications  of 
titanium  alloys.  A  complete  elastic-plastic  fracture  mechanics 
methodology  is  developed,  and  is  applied  to  example  models.  Topics 
where  further  research  is  needed  are  discussed  in  detail. \ 
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1.  INTRODUCTION 


1 .1  Background 

The  assurance  of  safe  and  reliable  structural  performance  of 
critical  components,  structures,  and  equipment  subjected  to  adverse 
loading  conditions  has  always  been  a  matter  of  vital  concern  to  both  the 
U.S.  Navy  and  the  Westinghouse  Electric  Corporation.  The  capability  to 
conduct  appropriate  structural  integrity  analyses  takes  on  an  added 
importance  when  new  equipment,  designs,  materials,  inspections  and 
fabrication  procedures  are  concerned. 

In  these  situations  there  is  little  or  no  service  experience  to 
rely  upon;  hence,  a  thorough  structural  integrity  analysis,  incorporating 
all  of  the  interacting  factors  must  be  included  as  a  major  element  in 
the  overall  plan.  Such  analyses  should  take  advantage  of  the  most 
advanced  technology  areas  that  are  applicable  to  the  situation  of 
concern;  in  this  case  modern  fracture  mechanics  technology  offers  a 
unique  and  directly  applicable  capability. 

Early  developments  of  fracture  mechanics  focused  on  plane-strain 
or  essentially  linear  elastic  fracture  conditions  (LEFM)  and  on  rela¬ 
tively  high  strength  brittle  materials  such  as  aircraft  structures, 
missile  cases,  gun  tubes,  etc.  Soon  the  technology  was  extended  to 
include  fatigue  and  stress  corrosion  crack  propagation.  Later  on, 
because  of  the  recognized  limitations  in  the  applicability  of  LEFM, 
effort  was  devoted  to  extend  the  fracture  mechanics  technology  to 
encompass  situations  involving  considerably  more  plasticity  than  is 
permissible  under  LEFM  conditions.  As  a  result  the  break  through  came 
in  the  form  of  the  path  independent  J-integral,  a  field  parameter 
analogous  to  K  in  LEFM.  The  general  usefulness  of  the  technology  has 
thus  been  extended  to  a  much  broader  spectrum  of  applications  and 


1 


‘  V *•  ** o'*  •*.  **.  >*.  •*.  1 


.v • 


.V  A 


'•/N  v; 


materials:  lower  strength,  higher  localized  stress  regions,  low  cycle 

fatigue  and  creep  controlled  crack  growth.  Even  more  recently,  the 
technology  has  taken  another  major  step  forward  with  the  advent  of  J 
resistance  curves,  tearing  modulus  concepts  and  tearing  instability 
models.  These  recent  developments  offer  the  capability  to  predict  the 
permisable  amounts  of  stable  crack  growth  in  the  ductile  temperature 
regime,  and  the  eventual  instability  conditions  for  the  catastrophic 
failure  of  the  structure  by  ductile  tearing  under  fully  plastic 
conditions.  More  importantly  these  recent  advances  in  technology  offer 
the  promise  of  enabling  the  design  of  structures  and  selection  of 
materials  so  as  to  avoid  any  possibility  of  failure  due  to  ductile 
tearing  instability. 

In  short,  fracture  mechanics  provides  engineers  with  a  powerful 
new  tool  for  more  effective  design  pertaining  to  structural  reliability. 

1.2  The  Program 

1.2.1  Objectives 

In  this  spirit,  the  Westinghouse  Electric  Corporation  conducted 
the  present  program  with  the  major  objective  of  assisting  the  Navy  in 
developing  and  applying  advanced  fracture  mechanics  technology  to  ensure 
structural  integrity  in  critical  applications  of  titanium  alloys.  The 
specific  objectives  were 

1)  Development  of  methods  for  assessing  Structural  Reliability 

2)  Respond  to  Specific  Navy  Concerns 

3)  Recommendation  of  methods  for  implementing  structural 
reliability  procedures 

1.2.2  Approach 

In  order  to  accomplish  the  mentioned  objectives  it  was  decided 
to  use  the  following  approach:  the  program  was  divided  in  two  phases. 


Phase  I  Assimilation  of  pertinent  Information 

In  this  phase  a  comprehensive  review  of  the  currently  available 
information  and  data  needed  for  applying  advanced  fracture  mechanics 
technology  to  structural  reliability  analysis  as  well  as  the  concerns  of 
the  Navy  was  conducted.  This  review  included  such  pertinent  information 
areas  as  loading  conditions  and  stresses  involved,  material  properties, 
fabrication  procedures  and  nondestructive  inspection  capabilities. 

Phase  II  Structural  integrity  analysis 

Phase  II  included  the  development  of  methods  to  assess  struc¬ 
tural  reliability.  After  considering  different  candidates  the  best 
possible  choice  was  proposed  using  the  latest  concepts  in  fracture 
mechanics  methodology.  The  method  was  applied  to  different  models. 

1.2.3  Organization  of  the  Present  Report 

This  report  is  organized  as  follows.  Phase  I  and  Phase  II  are 
discussed  in  detail  in  Sections  2  and  3  respectively.  In  Section  4,  the 
results  of  applying  the  proposed  methodology  to  specific  models  are 
shown.  In  Section  5  a  general  discussion  of  the  program  and  recommenda¬ 
tions  for  future  work  are  made.  In  Section  6,  the  help  of  Navy 
personnel  is  acknowledged.  Finally,  References,  and  Tables  and  Figures 
can  be  found  in  Sections  7,  and  8  respectively. 


2.  PHASE  I  ASSIMILATION  OF  PERTINENT  INFORMATION 


I 

2.1  Survey 

An  extensive  survey  was  conducted  to  gather  information  and 

data.  This  survey  included  visits  by  the  Westinghouse  Program  Team  to 

different  Navy  facilities.  The  main  points  addressed  in  the  survey  were 

i  a)  Navy  concerns.  An  adequate  characterization  was  made  of  the 

different  areas  that  concern  Navy  personnel. 

b)  Available  information  and  data.  A  list  of  areas  where 
experimental  data  exist  or  is  being  generated  today  was 
made . 


2.2  Results 
2.2.1  Navy  Concerns 

The  first  phase  of  this  program  is  "Assimilation  of  pertinent 
information  and  data."  To  do  this  several  meetings  were  held  between 
members  of  the  Westinghouse  program  team  and  Navy  technical  personnel  at 
a  number  of  different  Navy  facilities.  Also  Westinghouse  participated 
in  a  ONR  Ti-100  Workshop.  As  a  result  of  these  meetings,  many  areas 
were  identified  as  having  significance  to  the  problem  of  structural 
integrity  analysis  for  Ti-100  and  questions  were  raised  regarding  these 
areas;  a  list  of  those  follows. 

1)  Toughness  -  How  much  toughness  is  enough?  How  should  toughness 
be  characterized  for  structural  analysis? 

2)  Dynamic  Loading  -  What  is  the  fracture  behavior  under  dynamic 
loading  rates?  Can  the  parameters  used  for  conventional 
fracture  analysis  also  be  used  under  dynamic  loading? 

3)  Fatigue  -  Can  fatigue  to  failure  be  analyzed?  What  is  the 
effect  of  zero  to  compression  loading  on  fatigue  crack  growth 
analysis? 


4)  Low  Cycle  Fatigue  -  What  is  the  effect  of  this  on  crack  initia¬ 
tion  and  growth? 

5)  Failure  Criteria  -  No  failure  criterion  is  presently  used.  Can 
one  be  identified  particularly  for  analyzing  fatigue  to 
failure? 

6)  Crack.  Growth  Under  Sustained  Load  -  This  phenomenon  has  been 
observed  in  the  form  of  subcritical  cracking  and  delayed 
failure.  Is  this  related  to  environmental  influences,  creep  or 
time  dependent  fracture  toughness  behavior? 

7)  Effect  of  Prestrain  on  Fracture  Behavior  -  How  is  fracture 
toughness  and  other  fracture  behavior  influenced  by  an  initial 
prestrain? 

8)  Effect  of  Residual  Stresses  in  Welds  -  How  can  these  be 
measured  and  how  can  they  be  incorporated  into  a  structural 
integrity  analysis? 

9)  Thickness  Effects  -  How  are  these  accounted  for  in  structural 
analysis? 

10)  Scale  Models  vs.  Real  Structures  -  How  well  do  scale  models 
predict  fracture  behavior  in  real  structures?  What  models  are 
the  most  appropriate  ones  to  analyze? 

11)  Shop  Fabrication  vs.  Field  Fabrication  -  Are  properties  in 
field  welded  structures  as  good  as  those  in  shop  welded 
structures? 

12)  Explosion  Bulge  and  Tear  Tests  -  What  significance  do  these 
have;  Can  they  be  analyzed  using  advanced  fracture  approaches? 

2.2.2  Data  available  or  work  in  progress 

In  these  meetings  the  present  availability  of  experimental  data 
and  the  work  in  progress  in  different  areas  were  discussed.  A  list  of 
those  follows: 

1)  Dynamic  vs  Static  J-R  Curves 

The  effect  of  loading  rate  in  the  J-R  curve. 


2)  J-R  Curve  and  Toughness 

Material  toughness  data  and  the  possible  effect  of  variat!  ■ 
on  those  values,  i.e.,  thickness,  pres t rain ,etc . 

3)  Fatigue 

Data  fatigue  crack  growth  for  different  conditions. 

4)  Stress  corrosion  cracking 

Different  aspects  regarding  the  stress  corrosion  cracking 

5)  Sustained  Load  Crack  Growth 

Crack  growth  under  a  fixed  load  level  in  different  conditions. 

6)  Creep 

Studies  and  experimental  work  in  creep  and  creep  crack  growth 

7)  Environmental  Effects 

The  effect  of  environment  on  all  of  the  above  mentioned  areas. 

2.2.3  Result  of  the  Survey 

The  specific  results  of  the  survey  are  listed  above,  however 
some  general  conclusions  can  be  drawn 

1)  Qualitative  vs  Quantitative 

The  concerns  and  questions  raised  on  the  several  issues  are  not 
quantitative  rather  they  are  of  a  qualitative  nature.  The  question 
asked  is  whether  some  parameter  would  at  all  affect  the  result,  rather 
than  precisely  how  much.  This  shows  that  many  concerns  are  of  a  basic 
nature. 

2)  Diversity  of  Concerns 

The  survey  revealed  that  at  different  Navy  locations  and  among 
people  of  different  technical  functions  there  was  frequently  a  differ¬ 
ence  in  interests,  concerns  and  priorities. 

3)  Points  Address  Experimental  Work 

Most  of  the  points  of  concern  can  only  be  answered  by  addition.)] 

experimental  work. 


4)  Points  of  Concern  and  Basic  Research 

Several  points  that  were  raised  are  subjects  presently  being 
addressed  in  basic  fracture  research.  In  fact,  it  is  not  just  a 
question  of  how  to  apply  known  concepts  to  titanium,  or  how  titanium 
behaves  under  certain  conditions  compared  to  other  materials  extensively 
tested.  Instead,  several  points  of  concern  are  still  an  open  question 
in  other  areas  where  fracture  mechanics  is  much  more  advanced  and  has 
been  extensively  tested. 


3.  PHASE  II  STRUCTURAL  INTEGRITY  ANALYSIS 

3.1  Introduction 

3.1.1  The  J  Integral 

Recognizing  that  the  applicability  of  LEFM  was  necessarily 
limited,  researchers  in  fracture  mechanics  endeavored  to  extend  the 
technology  to  encompass  situations  involving  considerably  more 
plasticity  than  is  characteristic  of  LEFM.  The  breakthrough  came  in  the 
form  of  the  path  independent  J-integral ,  [ 1 ]  a  field  parameter  analogous 
to  K  in  LEFM.  As  a  direct  extension  of  linear-elastic  behavior,  J  can 
be  regarded  as  the  strength  of  the  stress-strain  fields  near  the  crack 
tip  for  nonlinear  elastic  materials. [2, 3, 4, 5]  Figure  1  illustrates  the 
crack  tip  stress-strain  field  with  J  as  the  single  parameter  characterizing 
the  strength  of  that  field. 

The  ductile  fracture  can  be  divided  into  separate  steps,  Figure 
2a.  Each  step  can  be  represented  on  a  plot  of  J  versus  ductile  crack 
extension,  Aa,  which  is  labeled  the  R  curve,  Figure  2b.  The  point  on 
the  R  curve  where  a  blunted  crack  tip  begins  to  tear  in  a  stable  manner 
has  been  labeled  Jjc»  Jjc  marks  the  initiation  of  the  ductile  cracking 
behavior  in  a  material  and  has  been  labeled  as  a  material  property.  Jjc 
is  limited  in  use  for  structural  fracture  evaluation  in  that  it  often 
provides  a  too  conservative  evaluation  of  toughness.  The  stable  crack 
advance  shown  in  the  R  curve  of  Figure  2b,  often  provides  useful 
structural  life  well  beyond  the  JIc  point.  The  R  curve  then  becomes  the 
representation  of  fracture  toughness  most  useful  for  evaluating 
structural  behavior. 
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3-1.2  R  Curve 


The  concept  of  a  crack  growth  resistance  curve  —  the  R  curve  — 
provides  a  useful  general  framework  for  understanding  the  relationship 
between  specimen  geometry  effects  and  material  behavior.  This  concept 
may  be  explained  in  the  following  manner:  Consider  a  cracked  body  being 
monotonically  loaded.  Any  fracture  mechanics  parameter  —  such  as  G,  J, 
K,  crack  opening  displacement  (COD),  or  crack  opening  angle  (COA)  —  can 
be  chosen  to  characterize  the  generalized  crack  driving  force.  The  crack 
may  start  to  grow  in  a  slow-stable  manner  above  a  certain  applied  crack 
drive  level.  The  R-curve  represents  a  locus  of  equilibrium  conditions 
where  the  crack  will  remain  stable  if  the  loading  is  stopped,  and  where 
the  applied  crack  driving  force  is  equal  to  the  crack  growth  resistance 
of  the  material;  i.e.,  Gj  or  J  =  R(MN-m/m^)  or  =  KR(MPa/m^^). 
Unstable  crack  growth  develops  when  the  crack  driving  force  increases  at 
a  greater  rate  than  the  crack  growth  resistance. 

The  R-curve  concept  has  been  successfully  developed  for  use  in 
the  linear-elastic  range,  and  the  characterizing  crack  drive  can  be 
either  or  G^.  Experiments  have  shown  that  in  the  dominant 
elastic/small-scale-plastic  range,  a  crack  growth  resistance  curve  can 
be  considered  a  material  property.  In  stating  this,  the  R-curve  is 
proven  to  be  independent  of  specimen  geometry  effects  such  as  initial 
crack  size,  specimen  width,  or  specimen  type.  Observed  effects  are 
dependent  on  the  material  thickness  and  test  temperature,  since  these 
variables  affect  the  basic  material  toughness  behavior. 

The  first  step  in  considering  elastic-plastic  and  fully  plastic 
R-curves  is  the  choice  of  a  suitable  parameter  to  characterize  the 
generalized  crack  driving  force.  J-integral  and  COD  concepts  were 
designed  specifically  for  use  up  to  the  crack  growth  initiation  point. 
After  significant  slow  stable  crack  growth  has  developed,  the  use  of 
these  elastic-plastic  methods  must  be  justified.  Hutchinson  and 
Paris [6]  showed  that  J  can  be  used  to  characterize  the  crack  growth 
process  if  the  remaining  ligament,  b,  in  a  J  test  specimen  is  large 


enough  so  that  a  region  of  proportional  strain  field  easily  encompasses 
the  local  crack  tip  nonproportional  strain  field.  They  defined  a 
material  based  length  parameter,  D,  which  they  expressed  as  a  function 
of  the  R-curve  slope  and  J  itself;  that  is,  D  =  J/(dJ/da).  And  then  the 
size  requirement  is  b  »  D  or  u  =  b/D  »  1.  The  required  size  of  the 
proportional  region  before  J  is  affected  has  not  yet  been  well  fixed. 

The  R  curve  based  on  J  as  the  characterizing  parameter  is 
justified  through  the  analysis  of  Hutchinson  and  Paris.  It  can 
therefore  be  used  to  assess  the  fracture  stability  of  structural 
components;  this  analysis  of  structural  stability  is  discussed  in  the 
following  section. 

3.2  Stability  Analysis 

Although  the  J-R  curve  concept  is  growing  in  acceptance  and  its 
use  becoming  more  common,  the  question  of  whether  a  crack  will  grow 
stable  or  unstable  remained  unresolved  until  very  recently.  It  was 
recognized  that  instability  results  from  a  lack  of  balance  between  the 
rate  of  increase  of  the  applied  drive  force  and  that  one  of  the  material 
resistance  to  crack  growth. 

Recently,  the  basic  implications  of  this  concept  were  further 
explored  by  Paris  et  al.[7-8]  and,  as  a  result,  it  was  demonstrated  that 
the  overall  characteristics  of  the  structure  play  a  major  role  in 
instability  and  its  effects  have  to  be  included  in  the  rate  balance 
mentioned.  In  this  work,  they  introduced  a  non-dimensional  quantity 
called  the  tearing  modulus,  T,  that  in  general  has  the  form: 

t  =  E  dJ 

T - 1  da  (1) 

a 

o 

where  E  is  the  elastic  modulus  and  aQ  is  the  flow  stress.  If  Equation 
(1)  is  evaluated  using  the  J-Resistance  curve  of  the  material,  the 
resulting  T  is  the  material  tearing  modulus  Tmat.  If  instead,  dJ/da  in 
Equation  (1)  is  calculated  as  the  rate  of  change  of  crack  drive  or  the 


applied  J,  per  unit  virtual  crack  extension,  with  the  condition  of  total 
displacement,  5tot,  kept  constant  (or  other  similar  conditions 
specified),  the  resulting  T  is  the  applied  tearing  modulus  Tapp.  And  so 
following  Ref erences 1 ,  instability  will  occur  when: 

^app  —  ^mat  (2) 

Using  the  condition  of  total  displacement  constant,  the  compliance  of 
the  structure,  C^,  is  introduced  into  the  analysis  the  Tgpp  becomes  a 
function  of  CM*  Consequently,  according  to  the  theory,  Equation  (2), 
instability  is  predicted  provided  the  values  of  Tmat  and  expressions  for 
Tapp  are  known.  In  their  original  work,  Paris  et  al.  developed  formulae 
for  Tapp  for  different  configurations  assuming  perfect  plasticity  and 
that  the  crack  grew  under  limit  load  conditions.  They  also  performed 
the  first  experimental  evaluation  of  the  theory.  In  tests  of  three 
point  bend  specimens  loaded  in  series  with  a  spring  bar  of  adjustable 
length,  the  compliance  of  the  system,  C^,  was  varied  from  test  to  test, 
producing  stable  or  unstable  behavior  in  complete  agreement  with  the 
theory. 

The  results  were  originally  shown  in  a  Tmat  vs  Tapp  plot. 

Figure  3.  Each  test  is  represented  by  its  Tmat  and  Tapp  values.  Open 
points  correspond  to  unstable  crack  growth  while  filled  points 
correspond  to  stable  crack  growth.  Ideally,  according  to  the  theory, 
all  open  points  (unstable)  should  lie  to  the  right  of  the  45°  line 
(region  where  Tflpp  >  Tmat),  while  the  filled  points  (stable)  should  lie 
at  the  left  of  the  45°  line  (region  where  Tmflt  >  Tapp).  The  excellent 
agreement  shown  in  Figure  3  was  the  first  experimental  verification  of 
the  tearing  instability  theory. 

3.2.1  Specimen  in  Series  with  a  Spring  Displacement  Controlled  Conditions 

Since  this  first  work,  significant  effort  was  devoted  to  extend 
the  same  concepts  in  different  directions.  To  better  illustrate  these 
concepts  consider  the  example  shown  in  Figure  4.  A  specimen  is  loaded 


in  series  with  a  linear  spring  of  constant  ,  simulating  tin- 

structure,  in  a  displacement  controlled  test. 

Hutchinson  and  Paris  developed  a  general  expression  for  T  tor 
this  configuration 
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where  6  is  the  displacement 

due  to  the 
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and  Cm  1s  che 

compliance  of 

the  spring.  All  the  terms  appearing  in  the  above  equation  are  calibra¬ 
tion  functions,  i.e.,  they  don't  bear  any  information  regarding  the 

material  response  to  crack  growth.  These  functions  can  be  obtained  from 

finite  element  analysis  or  experimentally  from  blunt  notch  specimen 

tests,  and  no  "real"  crack  growth  test  is  needed  for  their  determina¬ 

tion.  This  scheme  has  been  used[7-8j  to  obtain  Tapp  for  different 
configurations  of  practical  interest  and,  as  was  mentioned,  instability 
can  be  predicted  by  comparing  the  value  of  Tapp  obtained  from  Equation  (2) 
with  the  experimentally  obtained  Tmat. 

Recently  Ernst  et  al.[9,10,ll]  developed  general  formulae  to 
evaluate  both  Tmat  and  Tapp  directly  from  a  load-displacement  (P-6)  test 
record  giving,  for  the  above  case,  respectively 
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This  formulation  allows  a  direct  comparison  between  the  TmaC  vs  Tapp 
values.  The  instability  condition  can  now  be  found  directly  defining 

^T  =  Traat  -  Tapp  (6) 
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Instability  will  result  when 


AT  <  0 


or  using  Equations  (4)  and  (5). 
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AT  <  0  if  and  only  if 


-dP/d6  >  K  for  the  instability  condition 
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Thus  according  to  the  theory,  instability  will  ensure  when  At  <  0  and 
that  will  happen  if  and  only  if  -dP/d6  >  K^. 


a 


3.2.2  Alternate  Physical  Interpretation  of  the  Instability  Condition 

The  condition  for  instability  can  be  obtained  using  a  different 
approach. 

Consider  the  P-6  record  of  a  bend  specimen  tested  under 
displacement  control,  and  the  correponding  calibration  (nongrowing 
crack)  curves  as  shown  in  Figure  5a. 

It  is  assumed  that  the  fracture  process  is  described  in  terms  of 
the  J-R  curve;  that  is,  every  point  in  the  P-6  record  has  associated  a 
value  of  J,  J ^  and  a  value  of  a,  a^  which  are  connected  according  to  the 
J-R  curve. 

Suppose  now  that  an  identical  specimen  (same  a/W)  is  tested, 
this  time  in  series  with  a  spring  (as  described  before),  Figure  5b.  It 
can  be  seen  that  the  effect  of  the  spring  on  the  calibration  functions 
is  just  to  shift  every  point  in  Figure  5a  to  the  right  by  an  amount 
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where  PA  and  Pg  are  the  loads  at  Points  [A-A'J  and  [ B— B * ] ,  respec¬ 
tively.  Note  that  corresponding  points,  [A-A'],  [B-B' ] ,  etc.,  have  the 
same  value  of  J(J  depending  on  a/W  and  displacement  only  due  to  the 
crack,  6).  Thus  the  resulting  test  record  is  expected  to  go  through 
these  corresponding  points  [A*],(B'],  etc.  in  order  to  follow  the  J-R 
curve  as  before.  Combining  Equations  (10)  and  (11)  gives 


6B'  6A'  6B  6  A  +  CM(PB  “  PA^  (12 

is  important  to  note  that  for  the  portion  of  the  test  record  where  Pg  - 
PA  <  0  (dropping  part),  the  relative  distance  of  subsequent  points  is 
diminished  by  the  addition  of  the  spring 


6b  - 


5a  >  V 


-  V  if  PB  <  PA 


(13) 


In  fact,  if  enough  compliance  is  added,  this  relative  distance  can  even 
turn  out  to  be  negative;  that  is,  Point  [B'j  lying  to  the  left  of 
[A' ]  •  If  this  is  the  case  the  test  record  would  have  to  go  backwards 
(in  6)  to  pass  through  { B  * ]  In  order  to  follow  the  J-R  curve.  But  this 
Is  not  compatible  with  the  boundary  condition,  which  asks  for  a 
monotonically  increasing  displacement.  Thus,  the  test  record  gets  as 
near  to  Point  [B1]  as  it  is  allowed  to  (vertical  drop),  corresponding  to 
unstable  growth. 

The  conditions  for  stability  can  be  then  expressed  as 


6gt  -  6^1  <  0  unstable 


6 „ t  -  6 . ,  >  0  stable 


(14) 


Replacing  Equations  (10)  and  (11)  in  the  preceding  expression  gives 


CM  1  =  ^  <  -  dP/dS  for  instability 

This  condition  was  found  on  a  completely  general  basis  with  no  need  to 
mention  J  or  even  crack  length.  It  is  exactly  the  same  as  that  one  of 
Equation  (9)  found  from  tearing  instability  theory.  As  a  result  the 
condition  AT  <  0  is  a  necessary  and  sufficient  one  for  instability. 

Thus,  the  tearing  instability  theory  has  been  proven  to  be  always  valid. 


3.2.3  Specimen  in  Parallel  with  Spring  Load-Controlled  Conditions 

It  is  well  recognized  that  the  structural  member  is  normally 
constrained  by  the  structure  in  such  a  way,  that  neither  an  exclusively 
load-controlled,  nor  a  displacement-controlled  situation  is  realistic, 
but  more  likely  a  mixed  one  where  the  stiffness  of  the  rest  of  the 
structure  KM  has  to  be  taken  into  account.  In  the  last  section  it  was 
shown  that  by  adding  a  spring  in  series  with  the  specimen  in  a 
displacement  controlled  test,  the  structure  can  be  simulated  and  in 
fact,  the  whole  range  from  load  control  to  displacement  control  can  be 
covered  by  just  changing  the  stiffness  of  the  spring  %•  Nevertheless, 
that  is  not  the  only  possible  way  of  ranging  from  one  end  to  the 
other.  In  fact,  there  are  situations  where  the  structure  is  in  an 
intermediate  loading  condition  which  cannot  be  represented  by  the  above 
mentioned  model.  Namely  redundant  structures  where  several  members  are 
sharing  the  overall  applied  load:  stiffened  plates  or  pipes,  set  of 
structural  supporters,  cables,  etc.  The  stability  of  crack  growth 
problem  for  this  type  of  structures  was  the  subject  of  a  recent  paper  by 
Paris  et  al [ 12 ] ,  where  they  developed  the  so-called  Fracture  Proof 
Design  Concept.  According  to  this  concept  structural  parts  can  be 
designed  in  such  a  way  as  to  guarantee  stable  crack  growth  even  under 


load-controlled  conditions.  To  better  illustrate  these  concepts 
consider  a  specimen  in  parallel  with  a  spring  of  constant  K\j  simulating 
the  structure,  subjected  to  load  controlled  conditions,  as  shown  in 


Figure  6. 
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This  equation  is  the  general  expression  for  Tapp  for  a  specimen  in 
parallel  with  a  spring  (structure)  under  load  controlled  conditions .The 
instability  condition  can  be  obtained  by  comparing  Equations  (16)  and 
(17)  giving 


AT 


T  <  0 
app 


(18) 


if  and  only 
-dP/d6  >  Kj^ 
for  instability 


Following  a  similar  line  to  the  last  section,  the  instability  condition 
can  be  obtained  from  the  P-6  record  itself. 

Consider  the  P-6  record  of  a  specimen  tested  under  displacement 
control  as  shown  in  Figure  7a.  Suppose  now  that  an  identical  specimen 
is  tested  in  parallel  with  a  spring  of  compliance  CM  (0r  stiffness  %  = 
^).  The  displacement  underwent  by  the  specimen,  and  that  one  of  the 
spring  6^  are  same 


(19) 
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On  the  other  hand,  according  to  the  principle  of  equilibrium,  the  load 
held  by  the  specimen,  P,  plus  the  load  held  by  the  spring  PM,  give  the 
total  externally  applied  displacement  Ptot 

P  to  t  =  P  +  PM  (20a) 

As  a  consequence,  the  second  specimen  test  record  (Ptot>  $)  shown  in 
Fig.  7b,  will  be  just  the  result  of  shifting  every  point  of  the  original 
P-<5  record  up  in  load  by  an  amount 


PM  %  6  (20b) 

As  before,  if  two  generic  points  A  and  B  in  Fig.  7a  are  considered,  with 
coordinates  (6A,  PA)  and  (Sg,  Pg)  respectively,  the  addition  of  the 
spring  will  shift  them  to  (6A,  PA«)  and  (5g,  Pgt),  with 


PA'  "  PA  +  %  6  A 
PB'  =  PB  +  KM  5B 


(21) 


Combining  the  above  equations  one  gets  Equation  (22)  here 


(PB'  "  PA')  =  (PB  "  PA^  +  KM  (6B  '  «A>  (22 

The  above  expression  can  be  regarded  as  the  fundamental 
instability  equation  for  load  controlled  systems.  If  points  A  and  B 
were  in  the  raising  portion  of  the  P-5  record  (PA  >  Pg)  the  addition  of 
the  spring  causes  an  increase  in  the  load  difference  due  to  the  fact 
that  the  second  term  in  the  right  hand  side  of  Equation  (22)  is 
positive . 


<p„,  -  pAo  >  (pA  -  pB) 


> 


(23) 


If  the  points  A  and  B  were  in  the  dropping  portion  of  the  P-6  record 
(PA  >  Pg),  the  addition  will  also  cause  an  increase  in  the  load 
difference  for  the  same  reason  as  before. 

(PB.  -  PA.)  >  (PB  -  PA)  (24) 

In  fact,  if  the  second  term  in  the  right  hand  side  of  Equation 
(22)  is  big  enough  the  relative  difference  in  load  can  turn  to  be 
positive  Pgi  -  PA«  >  0,  and  that  means  an  always  increasing  curve  in 
load .  The  required  value  of  to  produce  this  behavior  can  be  obtained 
from  Equation  (22) 

(PB,  -  PA.)  =  (Pg  -  PA)  +  ^  («B  "  V  >  0 

P  -  P 
_  B  *A 

5.  -"5"  ^  Si  (25) 

B  A 

-dP/d6  <  1^ 

to  produce  a  monotonically  increasing  curve. 

Note  that  if  the  resulting  Ptot  -  6  record  is  a  monotonically 
increasing  curve,  it  is  irrelevant  whether  the  test  is  run  under  load  or 
displacement  controlled  conditions  i.e.,  the  test  will  be  stable 
regardless  of  these  conditions.  Thus,  Equation  (25)  represents  the 
fundamental  condition  for  stable  behavior  under  load  controlled 
conditions,  and  the  minimum  value  of  that  satisfies  it,  is  the 
minimum  value  of  stiffness  required  to  produce  a  monotonically 
increasing  curve  and  thus  prevent  instability.  As  a  result  the 
instability  conditions  are 

-  dP/d6  <  stable 

-  dP/d6  >  unstable 


(26) 


Note  that  this  condition  is  exactly  the  same  as  that  one  of 
Equation  (18)  derived  for  a  completely  different  situation;  in  both 
cases  stable  behavior  will  be  guaranteed  if  and  only  if 

-dP/d6  <  Km  (27) 

As  in  the  last  section,  note  that  this  condition  is  completely  general 
and  agrees  with  that  one  obtained  from  the  tearing  instability  theory. 
This  condition  obtained  from  the  tearing  instability  theory  agrees  with 
Equation  (18)  obtained  on  a  completely  general  basis,  based  on  nothing 
else  but  the  principle  of  equilibrium.  As  a  result,  the  tearing 
instability  theory  has  proven  to  be  always  valid  for  load  controlled 
conditions  in  the  sense  that  instability  will  occur  if  and  only  if  AT 
<  0.  Note  that  as  before,  if  any  other  parameter  X  =  X(6,  a)  was  used 
instead  of  J,  the  resulting  X-based  tearing  moduli,  TmatX  and  TappX, 
will  provide  an  automatically  validated  tearing  instability  condition. 

Finally,  it  is  to  be  noted  that  Equations  (15)  and  (27)  are  the 
same  although  they  represent  completely  different  situations.  The 
former  was  obtained  for  a  specimen  in  series  with  the  spring  under 
displacement  control,  while  the  latter  was  obtained  for  a  specimen  in 
parallel  with  the  spring  under  load  controlled  conditions. 

3.2.4  J-T  Diagrams 

In  the  last  sections  attention  was  devoted  to  the  development  of 
formula  for  Tapp  for  different  situations.  Although  since  the  original 
work  of  Paris  et  al.[7-8],  significant  effort  was  also  devoted  to 
experimentally  verify  the  tearing  instability  theory.  Joyce  et  al.[13], 
Vassilaros  et  al.[14]  and  Kamimen  et  al.[15]  among  others  conducted 
several  experimental  programs  to  do  this.  The  basic  philosophy  has  been 
the  same  as  that  one  of  the  initial  work  of  Paris  et  al.[7-8];  select  a 
certain  material-specimen  geometry  combination  and  run  a  series  of 
displacement  controlled  tests  of  specimens  in  series  with  a  spring  whose 
compliance  could  be  changed  from  test  to  test.  Values  of  Tmat  were 
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obtained  as  the  normalized  slope  of  the  J-R  curves  and  Tapp  values  were 
obtained  using  some  of  the  different  schemes [6-1 1 , 16 ] . 

As  was  mentioned  earlier  at  the  beginning,  the  results  were 

reported  in  a  Tfflat  vs  T app  plot  similar  to  that  of  Figure  3,  but  later 

on,  it  was  realized  that  in  general  neither  Tnat  nor  Tapp  were  constant 

ir.  a  given  test.  As  a  consequence,  the  results  were  shown  in  a  J-T  plot 

as  shown  in  Figure  8.  Slopes  of  the  J-R  curves  were  obtained  at 

different  J  levels  and  after  normalized  replotted  in  a  J  vs  Tma(; 

frame.  At  the  same  time,  using  some  of  the  schemes  mentioned  above 

values  of  were  obtained  for  different  J  levels  and  also  plotted  in 

app 

the  J-T  frame.  Every  time  that  the  two  curves  intercepted,  instability 
was  expected  to  occur  at  the  intersection  of  the  two  curves;  if  the  test 
were  stable,  the  two  curves  should  not  intersect.  If  that  was  the  case, 
instability  was  properly  predicted  by  the  tearing  instability  theory  and 
thus  the  latter  was  said  to  be  validated.  The  results,  as  predicted  by 
the  last  section  analyses,  have  proven  to  always  validate  the  tearing 
instability  theory. 

3.3  Methodology  for  Structural  Analysis 

In  the  previous  section  the  solution  to  the  ductile  fracture 
stability  problem  in  a  structure  was  given  with  the  model  of  a  test 
specimen  in  series  or  parallel  with  a  spring.  The  test  specimen 
represents  a  cracked  component  in  a  structure.  It  provides  the 
resistance  to  increased  crack  extension.  The  spring  in  series  repre¬ 
sents  the  driving  force  tending  to  pull  the  structure  to  instability  in 
a  displacement  limited  structure.  The  spring  in  parallel  represented 
redundant  components  resisting  instability  in  a  structure  under  applied 
loading.  The  problem  of  predicting  instability  in  the  structure  is 
simple  if  the  cracked  component  is  identical  in  geometry  to  the  labora¬ 
tory  specimen  used  to  generate  fracture  resistant  R  curves.  The  problem 
to  be  solved  in  providing  a  methodology  for  structural  analysis  is  that 
of  predicting  stability  when  the  cracked  component  represents  a  dif¬ 
ferent  geometry  from  that  of  the  test  specimen.  This  problem  can  be 


solved  without  conducting  a  test  for  each  new  geometry  if  two  sets  of 
information  can  be  generated  for  that  new  cracked  geometry.  They  are 
the  calibration  functions  for  the  geometry  which  relate  load,  dis¬ 
placement  and  J  as  a  function  of  crack  growth.  The  second  piece  of 
information  is  the  materials  resistance  to  crack  extension  which  is 
given  by  the  R  curve  for  the  cracked  structural  geometry. 

3.3.1  Calibration  Curves 

The  calibration  functions  can  be  obtained  from  the  load  versus 
displacement  behavior  for  the  non-growing  crack.  One  method  for 
obtaining  the  calibration  curves  for  a  given  cracked  geometry  is  by  the 
Shih-Hutchinson  estimation  procedures. [17-18]  This  procedure  allows 
displacement  and  J  to  be  calculated  as  a  function  of  load  and  crack 
length.  Inputs  needed  for  this  estimation  procedure  are  the  tensile 
properties  of  the  material  including  the  elastic  modulus,  E,  the  strain 
hardening  exponent,  n,  and  the  proportionality  constant,  a,  in  the 
Raraberg-Osgood  representation  of  hardening  behavior.  All  of  the 
constants  can  be  obtained  from  a  tensile  test  which  measures  the  true 
stress-strain  behavior. 

The  Shih-Hutchinson  estimation  procedure  has  been  applied  to  a 
number  of  different  cracked  geometries  using  the  finite  element 
calculation  method.  The  results  of  these  calculations  are  summarized  in 
a  Handbook  published  by  the  Electric  Power  Research  Institute. [ 19] 
Although  this  Handbook  does  not  cover  all  possible  geometries  of 
interest  and  the  accuracy  of  the  estimation  procedure  itself  has  not 
been  assessed  experimentally,  it  does  represent  a  viable  method  for 
obtaining  calibration  curves  for  many  different  types  of  cracked 
geometries . 

The  first  part  needed  to  apply  the  methodology,  that  of 
determining  calibration  curves  for  the  cracked  structure,  has  been 
solved  by  the  estimation  procedure.  The  results  are  available  in 
handbook  form  and  can  be  used  by  a  competent  engineer. 
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3.3.2  Material  Resistance  to  Crack  Growth 


The  second  piece  of  information  needed  to  apply  the  methodology 
is  the  R  curve  for  the  material  and  cracked  geometry  used  in  modelling 
the  structure.  This  is  a  more  difficult  piece  of  information  to 
obtain.  It  is  known  that  R  curves  may  be  geometry  dependent  so  that  one 
determined  from  a  laboratory  test  specimen  may  not  be  appropriate  the 
cracked  geometry  which  represents  the  structure.  What  is  needed  is  a 
way  to  represent  the  R  curve  as  truly  a  material  property  independent  of 
all  specimen  geometry  considerations.  One  step  toward  establishing  the 
R  curve  as  a  material  property  was  taken  in  specifying  limitations  for  J 
controlled  crack  growth. [6]  This  approach  stated  that  a  geometry 
independent  R  curve  can  be  developed  on  a  laboratory  test  specimen  which 
is  appropriate  for  the  cracked  structural  geometry  if  certain  limita¬ 
tions  can  be  met  for  both  the  laboratory  specimen  geometry  and  the 
cracked  structural  geometry. 

Unfortunately,  the  limitations  are  not  often  met  for  typical 
materials  used  in  building  structures.  A  significant  step  in  developing 
a  geometry  independent  R  curve  came  from  the  work  of  Ernst [20]  in  which 
he  suggested  a  modified  J  parameter,  J^,  which  could  be  used  to 
characterize  the  R  curve  behavior.  R  curves  plotted  with  JM  were  not 
subject  to  the  same  limitations  imposed  on  the  traditional  deformation 
J-R  curve  representation. 

Experimental  data  obtained  in  a  program  for  the  Electric  Power 
Research  Institute  were  used  to  show  the  adequacy  of  this  approach.  The 
program  is  described  in  detail  in  Ref erences [ 21-22 ] .  Basically,  it 
consisted  of  the  R-curve  testing  of  compact  tension  specimens  (CT)  of 
different  planar  sizes,  thickness  B  and  crack  aspect  ratio  a/W,  of  a 
A508  Class  2  steel  at  400°F. 

Out  of  the  large  test  matrix  some  specimens  were  selected  for 
this  study.  These  specimens  were  geometrically  similar.  All  of  them 
had  an  a/W  *  0.6  and  the  thickness  B  was  half  of  the  total  width  W  in  al 
cases.  The  overall  size  ranged  from  1/2T  (W  =  1  in)  to  10T  (W  =  20  in). 


The  resulting  J-R  curves  showed  consistency  at  the  early  stages  but,  as 
expected,  the  different  curves  start  deviating  when  the  crack  length 
increment  Aa  was  significant  compared  to  the  ligament  b. 


The  J-Aa  curves  were  then  converted  to  J-T  plots.  The  result  is 
shown  in  Figure  9.  Here  the  differences  between  specimens  are  more 
obvious  than  in  the  regular  J-R  curve  plot.  It  can  be  seen  for  lower 
Tmat  (higher  J)  the  "peeling  off"  of  the  smaller  specimen  curves 
occur.  In  fact,  as  was  mentioned  before  for  some  specimens,  Tmat  is 
close  to  zero  (even  negative  for  the  1/2T  ones).  The  same  data  were 
replotted  using  the  in  a  frame.  Using  for  JM  and  TM 


JM  =  JD  +  (1  +  0,76  b/W) 
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As  can  be  seen  in  Figure  10,  all  the  curves  collapse  into  one  when 
and  TMmar  are  used,  even  for  the  last  points  of  the  smallest 
specimens.  This  basically  means  that  JM  is  correlating  data  for 
situations  well  beyond  the  so-called  J  controlled  crack  growth  and  crack 
extension  up  to  30%  of  the  initial  remaining  ligament.  This  basically 
means  that  JM  is  a  better  parameter  to  use  for  describing  the  material 
resistance  to  significant  crack  growth,  and  that  its  use  for  R  curve 
representation  makes  maximum  meaningful  use  of  the  information  obtained 
from  small  specimens,  well  beyond  the  J  controlled  crack  growth  regime 
with  no  unnecessary  reduction  in  the  predicted  toughness  capability. 

3.3.3  Methodology 

When  the  two  components  needed  in  applying  the  methodology  are 
determined  (the  calibration  curves  and  the  R  curve  for  the  cracked 
geometry  representing  the  structure)  the  complete  response  of  that 
structure  to  an  applied  loading  can  be  determined.  From  this  response 


load  versus  displacement  plots  or  J  versus  T  diagrams  can  be  obtained. 
The  degree  of  stability  of  the  structure  and  load  capacity  can  be 
assessed . 

The  complete  methodology  is  represented  by  a  flow  diagram  in 
Figure  11.  The  structural  cracked  geometry  model  is  taken  with  the 
material  tensile  properties  to  determine  the  calibration  functions  for 
that  geometry.  The  R  curve  determined  from  a  laboratory  test  specimen 
is  represented  in  the  JM  versus  Aa  format.  This  represents  the  geometry 
independent  material  property  needed  to  assess  the  resistance  to  crack 
growth  in  the  structural  cracked  geometry.  The  calibration  curves  and 
the  R  curve  are  then  combined  to  produce  load  versus  displacement  or  J-T 
curves  for  the  model  geometry.  From  these  maximum  load  bearing  capacity 
and  stability  of  the  structure  are  assessed.  These  can  then  be  compared 
with  the  design  requirements  for  the  structure  to  determine  whether  the 
material  toughness  is  adequate.  Inadequate  toughness  would  require  that 
some  point  in  the  initial  input  be  changed.  This  could  be  a 
modification  of  design  requirements,  the  selection  of  a  material  with 
superior  toughness  or  even  a  reanalysis  of  the  structure  using  a  more 
refined  model  for  the  cracked  component  in  the  structure. 

3.4  Discussion 

The  methodology  presented  here  represents  in  principle  a 
complete  approach  to  assessing  the  load  bearing  capacity  and  stability 
characteristics  of  any  cracked  structural  geometry.  In  practice  there 
are  many  steps  in  the  methodology  which  need  further  assessment  and 
development . 

The  determination  of  the  calibration  functions  is  the  most 
advanced  part  of  the  method.  The  Handbook[19]  containing  the  Shih- 
Hutchinson  estimation  procedure  solutions  provides  the  information 
needed  to  develop  calibration  curves.  However,  the  number  of  solutions 
in  the  Handbook  is  presently  limited  and  may  not  exist  for  a  given 
structural  model.  Methods  for  adapting  a  solution  to  a  different 


geometry  or  combining  solutions  to  best  fit  a  structural  model  have  not 
been  developed.  Also  the  solutions  in  the  handbook  were  developed  by 
finite  element  computations  and  their  accuracy  has  not  been  verified 
independently  by  another  method  such  as  an  experimental  evaluation. 
Nevertheless,  the  Handbook  solutions  represent  the  best  approach 
presently  available  for  developing  calibration  curves  for  models  of 
cracked  structural  components. 

The  R  curve  evaluation  for  the  structural  model  represents  the 
most  difficult  part  of  applying  the  methodology.  The  method  developed 
by  Ernst [20]  representing  the  R  curve  with  the  modified  J,  JM,  is  the 
most  promising  approach  to  developing  a  geometry  independent  material 
property.  However,  this  approach  is  still  in  the  early  stages  of 
development  and  cannot  be  applied  with  complete  confidence.  Certainly 
this  approach  would  have  inherent  limitations  and  these  have  not  been 
determined. 

This  methodology  does  present  a  complete  approach  for  evaluating 
a  cracked  structure,  although  in  practice  there  are  some  steps  that  need 
further  work. 

First  of  all  it  is  to  be  noted  that  although  the  methodology  was 
presented  in  general,  its  use  was  exemplified  mainly  for  the  case  of 
slow  monotonic  loading,  assuming  that  the  material  properties  do  not 
change  with  time.  Obviously  there  are  other  mechanisms  of  deformation 
and  fracture  such  as  creep,  fatigue,  rate  dependent  monotonic  loading, 
environmentally  assisted  crack  growth,  and  all  the  possible  combinations 
of  the  above  mentioned.  Nevertheless,  the  methodology  presented  here, 
still  represents  a  complete  approach  for  evaluating  the  cracked  structure, 
if  the  correct  information  regarding  material  response  to  deformation 
and  crack  growth  for  the  particular  mechanism  is  input  in  the  central 
and  right  columns  of  Fig.  11  respectively.  In  other  words,  once  the 
mechanism  is  determined,  the  two  pieces  of  information  needed  are  still 
the  same:  the  calibration  functions  or  material  response  to  deformation 
and  resistance  curve  or  material  response  to  crack  growth,  obviously, 
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4.  APPLICATIONS 


4.1  Introduction 

In  the  last  section  the  elastic  plastic  fracture  mechanics 
methodology  was  presented.  This  scheme  represents  a  complete  approach 
for  evaluating  a  cracked  structure  in  terms  of  its  capacity  to  bear 
load,  displacement  and  energy  and  its  tendency  for  instability. 

It  was  mentioned  in  the  last  subsection  that  although  it  had 
been  somewhat  emphasi2ed  the  aspect  of  slow  monotonic  loading,  as  an 
example,  the  approach  can  be  used  for  any  other  mechanism,  provided  the 
necessary,  pertinent  information  regarding  material  response  to  both 
deformation  and  crack  growth  is  available.  When  more  than  one  mechanism 
is  present  a  study  regarding  how  the  different  effects  interact  between 
themselves  is  needed,  which  means  that  further  work  both  experimental 
and  analytical  has  to  be  done. 

One  of  the  outputs  of  the  presented  methodology  is  that  it 
clearly  shows  how  unnecessarily  conservative  the  design  concept  of  JIc 
really  is;  and  in  turn  allows  the  possibility  of  benefiting  from  the 
extra  toughness  of  the  material  developed  with  crack  growth.  This  is 
precisely  the  subject  of  this  section.  Basically  it  consists  of  a  very 
detailed  follow  up  of  the  methodology,  not  so  much  to  show  absolute 
numerical  values  but  rather  to  exemplify  with  a  sample  model  what  the 
method  can  do  in  terms  of  making  maximum  use  of  the  toughness  that  the 
material  really  has. 

Finally,  some  effort  was  also  directed  to  the  fatigue  crack 
growth  response  of  the  model. 
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4.2  Model 


As  was  mentioned  before,  the  first  step  in  the  methodology  is  to 
decide  upon  the  model  to  use.  This  choice  has  to  consider  two  aspects:  on 
one  hand  one  wants  the  model  to  represent  as  close  as  possible  the  real 
structure,  but  on  the  other  hand  there  are  usually  problems  regarding  the 
existence  of  solutions  (i.e.  calibration  functions  in  this  case),  the  full 
understanding  of  the  physics  of  the  problem  or  the  number  of  variables  to 
handle.  Thus  any  model  is  always  to  a  certain  degree  an  idealized 
representation  of  reality,  that  usually  neglects  part  of  it  due  to  the 
existence  of  the  mentioned  limitations.  In  this  case  the  structure  is  a 
cylinder  with  a  wall  thickness  W  =  4  in4^  and  inside  radius  R  =  200  in.  with 
circumferential  stiffners  at  a  distance  S  =  40  in  apart  subjected  to 
external  pressure.  The  cracks,  located  at  different  locations,  are 
considered  part  through  the  wall  thickness  and  are  subjected  to  multiple 
loads  that  are  nonproportional  and  produce  mixed  mode  loading  (I  and  II). 

The  simplest  case,  and  probably  the  most  important  one  is  that  of  a 
part  through  the  wall  crack  located  in  the  circumferential  direction, 
midway  between  two  stiffners.  The  crack  aspect  ratio,  depth  to  length  is 
(2c/a)  =  10,  Fig.  12.  Even  for  this  situation,  there  is  no  available 
finite  element  solution. 

Thus,  since  the  development  of  such  solutions  was  not  part  of  this 
program,  it  was  decided  to  use  the  closest  solution  available  in  the  EPH. 
The  choice  was  a  simple  edge  bar  subjected  to  three  point  bending.  The 
span  was  S  =  40  in.  and  the  width  was  W  =  4  in. 

Note  that  for  a  part  through  crack  of  an  aspect  ratio  as  big  as  10, 
it  can  be  considered  as  infinite  (2c/a  -*•  °»)  for  practical  purposes 
regarding  the  crack  drive  force  at  the  vortex  (deepest  point).  Also  note 
that  since  R/ t  *  50,  the  curvature  of  the  shell  can  be  neglected  for  that 

crack. 

The  fact  that  the  stiffners  are  attached  to  the  structure  (that 
is  not  a  simple  supported  situation)  can  be  simulated  in  principle  by 
considering  springs  at  the  ends  of  the  bar. 


4.3  Monotonic  Loading 


4.3.1  Tensile  Material  Properties 

It  was  assumed  that  the  material  follows  a  Ramberg-  Osgood  law 


given  by 
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where  oq  and  eQ  are  the  yield  strength  and  yield  strain  respectively. 
The  values  of  these  constants,  provided  by  the  Navy,  are 
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4.3.2  Elastic  Plastic  Formulation.  Calibration  Curves 

Once  the  geometry  has  been  selected  and  the  tensile  material 
properties  are  available,  the  next  step  in  the  methodology  of  Fig.  11  is 
to  obtain  the  calibration  functions  through  an  elastic-plastc 
formulation.  Here  the  Elastic  Plastic  Handbook  (EPH)[19],  that  follows 
the  original  scheme  of  Hutchinson  and  Shih  [17,18],  was  used.  In  the 
EPH  the  solutions  are  presented  as  follows.  First  consider  a  material 
that  obeys  a  constitutive  law  of  the  pure  power  type 
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For  this  material,  the  expressions  for  J  and  load  point  displacement  6 
(called  A  in  the  EPH),  or  fully  plastic  solutions  are  given  as 


J  =  a  o  e  b  h, (a/W,n)(P/P  ) 
o  o  1  *  o 
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(33) 
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where  P  is  the  load  per  unit  thickness,  a  is  the  crack  length,  b  is  flu 
remaining  ligament,  and  PQ  is  the  limit  load  per  unit  thickness  for  .1 
perfectly  plastic  material  (n  =  »).  The  expression  for  PQ  for  plane 


strain  is 


P  =  0.728  a  b“/L 
o  o 


whereas  for  plane  stress  PQ  is 


P  =  0.536  a  b  /L 
o  o 


On  the  other  hand  if  the  material  is  linear  elastic, 
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the  elastic  J,  Jel.  and  the  elastic  6,  can  be  expressed  as  [23] 

J  4>2  F2U/U, 
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with  VL  =  (5.58  -  19.57  (a/W)  ►  36.82  (a/b)  -  34.94  (a/W) 
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Finally,  when  a  material  obeys  .1  Rambcrg-'sgood  law  Eq  the  EPil 

estimates  the  total  value  of  I  and  '  bv  adding  the  linear  elastic 
expression  plastically  connected  and  the  fully  plasic  solution.  The 
latter  consists  of  considering  an  effective  crack  length  which  is 


4.3.3  Material  Response  to  Crack  Growth 

The  material  response  to  crack  growth  was  characterized  by  the 
J-R  curves  provided  by  the  Navy,  for  the  three  material  conditions  0,  1 
and  3%  prestrain.  These  curves  were  obtained  using  IT  compact  specimens 
(W  =*  2”),  under  static  loading  at  room  temperature.  Moreover,  the 
curves  were  replotted  using  the  modified  J,  JM,  concept  to  give  JM  vs  Aa 
curves.  The  expression  for  Jj^  used  (valid  for  CT  specimens)  was 


(1  +  0.76  b/w)  J  . 
J  +  /*  - _ - — — 


where  Jpj  is  the  plastic  part  of  J  (total  minus  pure  elastic).  These 
plots  are  shown  in  Figures  13,  14  and  13. 

4.3.4  Computer  Program 

The  next  step  in  the  methodology  of  Fig.  11,  is  to  combine  the 
information  regarding  deformation  (central  column)  and  crack  growth 
(right  column)  in  the  computer  program  that  constructs  the  P-6  record  or 
J-T  diagrams  for  the  untested  geometry.  The  mechanics  of  the  program 
developed  for  this  purpose  was  as  follows. 
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First,  the  calibration  functions  given  by  the  F.PH,  the  -  Aa 
curve  and  the  expression  connecting  J  and  JM  for  the  selected  untested 
geometry  were  stored  into  the  program.  That  is 


J  =  J(P,a) 


6  =  6(P,a) 

-  R  curve 


Calibration  Curve 


pairs  (JM,  Aa) 


JM  =  J  +  J3  da 

M  J  a  b 

o 


Expression  for  J.,  for  3  point  bend 
specimen.  (42) 


After  some  value  of  the  initial  crack  length  was  selected,  the 
program  increases  the  load  P  from  0  to  the  point  where  the  corresponding 
J  matches  the  first  value  of  JM  in  the  R  curve  (Aa  =  0).  Then  the 
program  increments  aQ  by  Aa  and  with  this  new  crack  length  a  =  aQ  +  Aa, 
increases  the  load  P  from  0  to  the  point  where  the  corresponding  J  is 
such  that  produces  a  (through  Eq  (42)  that  matches  the  J^-R  curve  for 
that  Aa.  The  process  goes  on  until  the  whole  -  Aa  curve  is  covered. 

Every  time  that  the  correct  J  is  found  the  corresponding  values 
of  P  and  6  (or  J  -  Tapp  -  Toat,  or  JM  -  TMapp  -  TMmat)  are  saved  and 
reploted  to  build  the  P  -  6  record  of  the  untested  geometry  (or  J-T 
diagram),  which  follows  the  -  Aa  curve. 

4.3.5  Results  and  Discussion 

The  results  of  the  application  of  the  methodology  to  the  example 
model  are  shown  in  Figures  16  through  23  and  in  Table  1. 

The  first  three  figures,  16,  17  and  18  show  the  P-6  records  for 
aQ  ”  0.5;  1.0;  1.5  and  2.0  inches  using  the  J^-Aa  curve,  for  the  three 
material  conditions  0,  1  and  3%  respectively. 


The  three  subsequent  figures,  19,  20  and  21  show  the  same 
plotting,  using  this  time  the  deformation  J-Aa  curve. 

In  Figure  22,  the  P-A  records  for  aQ  =  1.5  in  and  the  different 

material  conditions  are  compared  using  the  -  Aa  resistance  curve. 

Figure  23  shows  plots  of  J  vs  Tapp  for  the  model  with  springs  in 

parallel.  The  value  of  the  spring  constants  per  unit  thickness  were 

Kj^j  =  10;  50  and  100  k  lbs/in2.  The  value  of  aQ  was  1.5  in.  A  schematic 

of  J  vs  T„„.  is  also  shown, 
mat 

In  Table  I,  the  results  for  all  the  cases  run  are  summarized  in 
terms  of  loads  P  and  values  of  J.  Three  values  of  the  load  are  compared: 
the  one  at  Jjc»  and  the  maximum  attained  ones  using  the  J-Aa  and  JM  -  Aa 
curves.  Note  that  the  maximum  attained  load  is  not  necessarily  the 
maximum  load  that  the  structure  can  withstand. 

Two  values  of  J  are  compared:  the  Jlc,  and  the  last  value  of  JM 
used  from  the  R  curve  i.e.  the  JM  corresponding  to  the  last  P  obtained 
through  the  computer  program. 

From  the  figures  it  can  readily  be  seen  that  in  most  cases,  the 
available  information  in  terms  of  material  resistance  to  crack  growth 
did  not  allow  the  program  to  attain  the  maximum  load,  Pmax»  that  the 
structure  can  hold.  As  a  consequence,  the  slope  of  the  decreasing 
portion  of  the  P-6  record  dP/d6,  needed  for  stability  calulations,  could 
not  be  obtained  either. 

This  basically  means  that  the  resistance  curves  provided  were 
too  short  (in  Aa)  to  obtain  points  in  the  P-6  record  beyond  Pmax-  1° 
other  words,  the  last  point  used  from  the  J  or  -  Aa  curves  was  not 
sufficient  to  get  to  Pmax*  Thus  the  values  of  J  or  Jj^  corresponding  to 
Pmax  are  even  bigger  than  the  last  ones  used,  and  those  are  already 
significantly  bigger  than  JIc-  This  can  be  appreciated  better  by 
considering  Table  I.  For  every  case  run,  the  load  at  JIc  is  obtained 
with  the  maximum  attained  value  by  using  the  -  Aa  curve.  It  can  be 
concluded  that  the  load  that  the  structure  can  hold  is  at  least  up  to 
20%  bigger  than  that  predicted  by  the  use  of  JIc’  On  the  other  hand,  in 


the  last  two  columns  of  the  table  the  value  of  Jj  is  compared  with  that 
of  for  the  maximum  attained  load.  It  can  be  concluded  that  the 
latter  are  bigger  by  a  factor  of  2  to  3.  This  fact  clearly  shows  how 
important  the  extra  toughness  developed  with  crack  growth  is.  The 
toughness  available  in  the  materials  is  really  much  bigger  than  . 

A  word  is  in  order  here,  regarding  the  R-eurves  of  Fig.  13,  14, 
15.  These  curves  were  obtained  using  IT  CT  specimens,  thus  the  amount 
of  crack  extension  is  roughly  10  to  15%  of  the  initial  uncracked 
ligament.  This  means  that  the  investigators  who  obtained  these  curves, 
did  the  right  thing  in  stopping  the  tests  and  those  points,  since  it  was 
believed  at  the  moment  that  the  amount  of  crack  extension  should  not 
exceed  10%  of  the  ligament  in  order  to  keep  the  crack  growth  under  the 
so-called  J  controlled  conditions. 

This  point  emphasizes  the  advantages  that  provides.  In  fact, 
as  shown,  JM  can  be  used  even  for  amounts  of  growth  of  30%  of  the 
initial  ligament.  This  means  that  much  longer  (in  Aa)  curves  can  be 
obtained  from  the  same  specimen  size.  Thus  allowing  the  evaluation  of  a 
much  larger  structural  parts  (like  in  this  case). 

On  the  other  hand,  by  looking  at  Fig.  22,  it  can  be  concluded 
that  the  material  condition  did  not  make  much  difference  in  terms  of  the 
obtained  P-6  records.  Although  it  is  emphasized  here  that  the  records 
shown  are  only  a  portion  of  the  complete  ones,  and  presumably  the 
differences  are  more  obvious  in  terms  of  load,  displacement,  and  energy 
when  longer  R-curves  are  used  and  thus  complete  P-6  records  are 
considered.  At  the  same  time,  as  can  be  seen  by  considering  Eq  (49),  ,J 
is  proportional  to  P^  in  the  elastic  regime  and  to  Pn+^  in  the  fully 
plastic  one.  Thus,  for  all  cases,  changes  in  J  represent  much  smaller 
changes  in  P.  Thus,  moderate  differences  in  resistance  curves  (due  to 
material  condition  for  example)  are  translated  into  much  smaller  ones  in 
P-6  records. 

Finally,  following  the  spirit  of  Section  3.2  the  model  was 
considered  loaded  in  parallel  with  springs  of  different  constants  to 


simulate  the  rest  of  the  structure.  The  intention  of  this  part  was  to 
perform  an  instability  analysis  following  Section  3.2,  by  measuring  the 
unloading  part  of  the  P-6  record,  dP/d6,  and  show  that  if  a  spring  of 
enough  stiffness  issued  in  parallel  with  the  model,  the  crack  growth  can 
be  stable  even  under  load  controlled  conditions,  past  maximum  load. 

Unfortunately,  as  was  extensively  discussed,  the  J-R  curves 
provided  where  not  long  enough  to  even  reach  Pmax,  thus  the  slope  of  the 
decreasing  part  of  the  P-6  records  could  not  be  measured.  In  terms  of 
J-T  diagrams  this  means  that  the  material  line  is  not  long  enough  to 
cross  the  applied  line.  In  any  event,  as  an  example,  J  vs  Tapp  values 
were  calculated  for  aQ  =  1.5  in  using  spring  stiffness  (per  unit 
thickness)  of  10,  50  and  100  k  lbs/in^,  a  schematic  material  line  (the 
real  on  does  not  cross  the  J-Tapp)  is  also  shown. 

In  summary,  it  can  be  concluded  that: 

1)  Available  J-R  curves  are  not  long  enough  (in  Aa)  to 
characterize  crack  growth  behavior  in  the  model.  (due  to  a  restriction 
on  Aa  to  have  J-controlled  crack  growth  conditions). 

2)  The  newly  developed  JM,  allows  one  to  greatly  extend  the 
range  of  valid  data.  As  a  result  its  use  is  recommended  to  make  maximum 
use  of  the  information  obtained  from  laboratory  specimens. 

3)  An  elastic  plastic  methodology  can  be  used  to  assess 
structural  integrity  in  a  model,  determining  its  characteristics  to  hold 
load;  displacement  and  energy  input  and  its  degree  of  instability.  4) 
The  presented  methodology  makes  full  use  of  the  real  toughness  of  the 
material  by  considering  also  the  extra  one  developed  with  crack 
growth.  Specifically  in  this  particular  model  used: 

a)  the  load  carrying  capacity  increases  at  least  5  to  30%, 

b)  the  toughness  increases  by  a  factor  of  at  least  2  to  3,  by 
using  the  present  methodology  as  opposed  to  the  unnecessarily 
conservative  J^c  design  criterion. 


4.4  Cyclic  Loading 


In  the  last  sections,  the  case  of  monotonic  loading  was  treated 
extensively.  This  usually  represents  the  case  of  accident  conditions 
where  loads  are  very  high  and  small  duration. 

On  the  other  hand,  also  important  is  the  problem  of  service 
conditions  where  loads  are  not  that  high  but  are  cyclically  applied. 

For  the  cyclic  loading  case,  elastic  plastic  fracture  mechanics 
methods  can  also  be  applied,  and  are  really  in  order  for  situations 
where  the  loads  are  high  enough  to  produce  extensive  yielding.  If  this 
is  the  case,  as  shown  by  Begley  and  Dowling[24],  the  parameter  to  use 
is  AJ,  which  in  turn  is  general  enough  to  cover  Linear  elastic 
situations  also. 

Thus,  neglecting  crack  closure  effects,  for  load  excursions  from 
zero  to  maximum  (R  =  0)  AJ  can  be  obtained  by  using  Eq  (40),  simply 
replacing  P  by  AP.  Here,  as  an  example,  this  methodology  was  applied  to 
the  model  detailed  above.  Basically,  a  computer  program  was  developed 
that  calculated  the  number  of  cycles  N  required  for  a  certain  initial 
crack  aQ  to  get  to  the  critical  a,  ac,  for  which  J  =  J^r  (at  P  x) . 

Thus  the  corresponding  material  properties,  J^c  and  the  expression  for 
da/dN  in  terms  of  J  provided  by  the  Navy  were  entered  into  the  program 

JT  =  580  lbs  in/in^ 

ic. 

da/dn  =  C  (AJ)r  ^3^ 

with 

C  =  io~7 
r  =  2 

when  a  is  measured  in  inches  and  J  in  lbs  in/in  . 

A.s  mentioned,  the  number  of  cycles  to  failure  were  calculated 
for  different  initial  aQ.  The  results  are  shown  in  Fig.  24.  As  1  final 
comment  it  can  be  said  that  the  material  response  to  fatigue,  Eq .  (43), 


5.  FINAL  DISCUSSION  AND  RECOMMENDED  RESEARCH 


In  the  previous  sections  the  elastic  plastic  fracture  mechanics 
to  assess  structural  reliability  of  cracked  structures  was  developed  In 
general  and  applied  to  specific  examples.  The  method  incorporates  new 
EPFM  concepts  through  the  use  of  a  J  based  R  curve  and  tearing  insta¬ 
bility  analysis.  The  central  feature  of  the  methodology  is  a  model  of 
the  cracked  structural  geometry  which  is  combined  with  a  spring  in 
series  or  parallel  for  the  stability  analysis.  The  method  is  based  on 
the  analysis  of  the  response  of  the  structure  to  deformation,  labeled 
the  calibration  curves,  and  the  response  of  the  material  and  structure 
to  crack  advance,  the  R  curve. 

The  calibration  curves  are  evaluated  through  the  Shih-Hutchinson 
estimation  procedure.  Results  for  this  procedure  are  available  in  an 
EPRI  published  Handbook.  The  R  curve  is  evaluated  as  a  geometry  inde¬ 
pendent  material  property  using  the  modified  JM  approach  developed  by 
Ernst.  By  combining  the  R  curve  and  the  calibration  curves  a  complete 
description  of  the  fracture  behavior  can  be  given  for  the  structural 
model.  Such  things  as  maximum  load  bearing  capability  and  structural 
stability  can  be  assessed. 

The  results  obtained  show  that  levels  of  toughness  between  2  and 
3  times  more  than  the  JIc  point  can  be  attained  before  reaching  the 
instability  point.  In  terms  of  load  carrying  capacity  the  benefit  of 
this  method  represents  at  least  an  increase  of  20% - 

The  methodology  presented  in  the  first  section  is  not  necessarily 
limited  to  slow  raonotonically  increasing  load  at  room  temperature.  In  fact 
it  can  be  used  when  other  mechanisms  of  deformation  and  crack  growth  are 
operating,  provided  the  correct  material  response  to  deformation  (central 
column  of  Fig.  11)  and  to  crack  growth  (right  column  of  Fig.  11)  are 
available . 
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In  principle,  situations  of  obvious  interest  for  the  Navy  as 
creep,  dynamic  loading,  Mode  II  loading,  among  others,  can  be  treated 
using  this  methodology.  However,  there  are  some  areas  of  uncertainty 
that  further  work,  both  analytical  and  experimental,  is  needed. 

The  Handbook[19]  containing  the  Shih-Hutchinson  estimation 
procedure  solutions  provides  the  information  needed  to  develop  calibr¬ 
ation  curves.  However,  the  number  of  solutions  in  the  Handbook  is 
presently  limited  and  may  not  exist  for  a  given  structural  model. 

Methods  for  adapting  a  solution  to  a  different  geometry  or  combining 
solutions  to  best  fit  a  structural  model  have  not  been  developed.  Also 
the  solutions  in  the  handbook  were  developed  by  finite  element  computa¬ 
tions  and  their  accuracy  has  not  been  verified  independently  by  another 
method  such  as  an  experimental  evaluation.  Nevertheless,  these  Handbook 
solutions  represent  the  best  approach  presently  available  for  developing 
calibration  curves  for  models  of  cracked  structural  components.  The  R 
curve  evaluation  for  the  structural  model  represents  the  most  difficult 
part  of  applying  the  methodology.  The  method  developed  by  Ernst [20] 
representing  the  R  curve  with  the  modified  J,  JM,  is  a  promising 
approach  to  developing  a  truly  geometry  independent  material  property. 
However,  this  approach  is  still  in  the  early  stages  of  development  and 
cannot  be  applied  with  complete  confidence.  Certainly  this  approach 
would  have  inherent  limitations  and  these  have  not  yet  been  determined. 
Other  geometries  should  be  tested,  and  those  limitations  should  be 
explored. 

The  effects  of  different  variables  affecting  the  J-R  curve  such 
as  loading  rate,  sustained  load,  cyclic  loading  (very  high  cyclos)  among 
others  and  the  combination  of  them  should  be  studied. 

In  conclusion,  the  methodology  developed  does  present  a  complete 
approach  for  evaluating  a  cracked  structure.  However,  there  are  certain 
areas  where  further  work  both  experimental  and  analytical  is  needed. 
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Fig.  2b-Four  steps  of  ductile  fracture  process  on  an  R  curve 
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A  508  Class  2A 
Temp  400°  F 
Compact  Specimens 
a/W  ~  0. 6 
B  =  Proportional  (W/2) 


Fig.  16-Load  per  unit  thickness  vs  total  deflection,  using  JM  for  different  a  ( 0%  prestrain) 
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Fig.  22-Load  per  unit  thickness  vs  total  deflection  using  for  aQ  =1.  5  and  different  prestrain 
conditions  J 
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Fig.  23— J -T  diagram  for  different  spring  compliances 


